Microscopic derivation of Frenkel excitons in second quantization 
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Starting from the microscopic hamiltonian describing free electrons in a periodic lattice, we derive 
the hamiltonian appropriate to Frenkel excitons. This is done through a grouping of terms different 
from the one leading to Wannier excitons. This grouping makes appearing the atomic states as a 
relevant basis to describe Frenkel excitons in the second quantization. Using them, we derive the 
Frenkel exciton creation operators as well as the commutators which rule these operators and which 
make the Frenkel excitons differing from elementary bosons. The main goal of the present paper is 
(3JT)' to provide the necessary grounds for future works on Frenkel exciton many-body effects, with the 

[""t , composite nature of these particles treated exactly through a procedure similar to the one we have 

. recently developed for Wannier excitons. 
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The absorption of photons in a dielectric solid can lead to delocalized excitations called excitons. These excitons 
are essentially of two types: the Wannier excitons^ and the Frenkel excitons^. 

Wannier excitons are formed in inorganic semiconductors. The relative motion of the electron and the hole from 
which the excitons are made, encompasses hundred of unit cells. This leads to a small binding energy (~ 10 meV) and 
a large Bohr radius (~ 10 nm). As a consequence, Wannier excitons start to interact at relatively low densities, giving 
rise to a large variety of many-body effects associated with optical nonlinearities, which makes these conventional 
Q semiconductors used in today's technologies. 

The second type of exciton, known as Frenkel exciton, is commonly found in organic crystals. These crystals are of 
potential importance for future electronic devices, which makes them under current intensive studies&^&i^. Frenkel 
excitons are formed with electron and hole localized on a small scale, of the order of a single molecular block (~ 1 
■ nm), the typical binding energy being of the order of 1 eV. Due to interactions between blocks, these single molecule 
I excitations are transferred from site to site, giving rise to a wave known as Frenkel exciton. 

These pictures can be qualitatively understood by noting that, in conventional (inorganic) semiconductor, the 
relative dielectric constant is rather large (~ 10), which makes the screening of the interaction between carriers quite 
strong. As a result, the attraction between electrons and holes is weak, which explains the large extension of their 
relative motion wave function. On the opposite, the small relative dielectric constant of molecular crystals (~ 1) leads 
to a strong electron-hole attraction which localizes the pair on a system unit cell. 

The separation between Wannier and Frenkel excitons, of course, is not very sharp. Recent works^ have shown 
the limitations of the simple Frenkel picture for excitons in organic semiconductors and the necessity to introduce 
"charge-transfer excitons" in which the electron and the hole are located in different sites, with similarity to Wannier 
excitons for which the distance between electron and hole is large compared to the ion-ion distance. Along the same 
line, recent progresses in device nanofabrication now allow one to obtain organic/inorganic semiconductor structure 
' in which the hybridization between Frenkel and Wannier excitons can be produced^. 

Due to their small exciton Bohr radius, interactions between Frenkel excitons are expected to occur at a much 
larger density than the one at which many-body effects between Wannier excitons start to be noticeable. This barely 
comes from the fact that the dimensionless parameter which controls these many-body effects is 

?7 = N(a x /L) D . (1.1) 

where a x is the exciton Bohr radius, L the sample size, D the space dimension, and N the exciton number - the 
exciton density being n = N/L D . 

Most likely, as for Wannier excitons, interactions between Frenkel excitons are going to be of importance in electronic 
devices constructed with organic semiconductors. This is why a correct treatment of these interactions is highly 
desirable. Being made of indistinguishable carriers, Frenkel excitons, like Wannier excitons, are not well defined 
objects, which makes the interactions between excitons not possible to identify properly. As a direct consequence, one 
cannot describe these interactions through a potential, as usually done by lack of a correct procedure. 
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Over the last few years, we have developed a new many-body theory 9,10 for Wannier excitons, in which the composite 
nature of the particles is treated exactly. We have shown that Wannier excitons predominantly interact through 
the Pauli exclusion principle which exists between their fermionic components. This exclusion gives rise to very 
many carrier exchanges between excitons which are nicely visualized through Shiva diagrams^ 1 -, rather different from 
Feynman diagrams due to the composite character of the particles. All our works on Wannier excitons end with 
the same conclusion: it is not possible to replace composite excitons by elementary excitons, as commonly done 
through sophisticated bosonization procedures 1 ^; either one misses terms as large as the ones kept, or, in optical 
nonlinear effects, one even misses the dominant terms 1 ^. This can be readily seen from a dimensional argument: the 
Pauli scatterings associated with carrier exchanges are dimensionless, while the Coulomb scatterings are energy-like 
quantities; so that they have to appear with an energy denominator, which can only be a photon detuning. This 
makes these Coulomb terms completely negligible in front of the pure exchange terms - missed with bosonized excitons 

- when unabsorbed photons have a large detuning. 

The development of a similar procedure for Frcnkcl excitons requires to settle a second quantization formalism for 
these excitons on a clean basis in order to possibly keep their composite boson nature exactly all over the calculations. 
In contrast, the works on interacting Frenkel excitons, we have up to now see n 14 ' 15 , contain a potential written as 

JmnBJnBn, where B\ creates an excitation on site n. Besides the fact that interactions between excitons cannot 
be written as a potential due to the composite nature of these particles, by writing such a potential in terms of B^ 
only, one obviously losses the composite electron-hole nature of this excitation which exists in the exact potentials 
between electrons and holes (as seen below from Eqs. (4.6,13). This is exactly this composite nature that we want to 
treat properly, since we know its importance in the case of Wannier excitons. 

In this first work on Frenkel excitons, we propose a microscopic approach to the description of these excitons based 
on a second quantization formalism, starting from the hamiltonian of free electrons in a periodic lattice. Through a 
grouping of terms different from the one leading to Wannier excitons, we introduce the atomic states as a physically 
relevant one-electron basis for Frenkel excitons and we rewrite the system hamiltonian in terms of electrons and 
holes localized on atomic sites in order to have a precise description of the interactions. Since we use a second 
quantization scheme in terms of electrons and holes separately - and not in terms of their product B^ n as usually done 

- our approach automatically takes into account the fermionic composite nature of the particles forming the Frenkel 
excitons. This is going to be of crucial importance for a proper study of many-body effects involving these excitons. 
The present preliminary work actually provides the necessary grounds for further works on Frenkel exciton systems. 
In a forthcoming publication, we are going to use this second quantization formalism to derive the Coulomb and Pauli 
scatterings of two Frenkel excitons: these are the elementary scatterings on which all many-body effects dealing with 
excitons are based. We will then use these scatterings to calculate the ground state energy of N Frenkel excitons in 
the low density limit, a physical quantity of basic relevance. 

The present paper is organized as follows. In Section II, we start with the first quantization description of a 
periodic system made of interacting ions and electrons. We discuss the conceptual difference between Wannier and 
Frenkel excitons which leads to a different grouping of terms in the hamiltonian. We introduce the atomic states as 
the appropriate one-electron basis for the problem when the tight-binding approximation, which neglect the overlaps 
between the atomic wave functions for electrons on different sites, is valid. We also discuss the conceptual difficulty 
associated to the atomic basis compared to the free electron and hole basis used in the case of Wannier excitons. In 
Section III, we derive the semiconductor hamiltonian appropriate to Frenkel excitons in second quantization using 
this atomic basis. In Section IV, we switch to holes and we reduce the hamiltonian to the terms which conserve the 
number of electron-hole pairs. We then discuss all these terms with a particular attention to the one responsible for 
the excitation transfer from site to site. In Section V, we identify the lowest excited states of the hamiltonian in the 
absence of interactions between sites and we show that they form a degenerate subspace. This degeneracy is split 
to give rise to Frenkel excitons by the intersite interactions. They are introduced in Section VI, which is devoted to 
the precise derivation of the creation operators for Frenkel exciton and the commutation rules which govern these 
operators. The precise handling of these commutators are at the basis of the many-body theory we are going to 
construct. We show that, like Wannier excitons, Frenkel excitons are composite bosons, their commutation rules 
differing from the ones of elementary bosons due to the fermionic nature of the electrons and holes forming these 
excitons. In Section VII, we collect the main results of this paper and conclude. 

II. FIRST QUANTIZATION DESCRIPTION 

A. Semiconductor Hamiltonian in first quantization 

Let us consider N s electrons with charge — |e|, located at and N s ions with charge + |e|, located at R„, with i 
and n running from 1 to the number of sites N s . In this first work on Frenkel excitons, we are going to forget all 
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spin degrees of freedom for the sake of simplicity. This physically corresponds to have all the electrons with the same 
spin. This also means that we drop all degeneracies coming from the orbital part of the electronic levels. These spin 
and orbital degrees of freedom generate very interesting polarization effects. They, however, lead to heavy notations 
which are wise to avoid in a first work. 

The semiconductor Hamiltonian can be written as 

H — Hkin ~t~ Ve—ion ~l~ Vee ~t~ Vion — ion- (^--0 

The one-body operator Hkin describes the electron kinetic energy 
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where m is the free electron mass. The second term of Eq. (2.1), which describes the electron-ion Coulomb attraction, 
also is a one-electron operator. It reads 

N s N s 2 

i=i „=i i r » "-"i 

The third operator V ee describes the Coulomb repulsion between electrons. This two-body operator is given by 

"-^E^- (") 

The last term, Vi on -i on , which describes the Coulomb interaction between ions, is a constant with respect to the 
electron motion. It is however necessary to keep it in order to work with the hamiltonian of a fully neutral system. 
This is required for the convergence of the Coulomb terms in the large sample limit. We will see below the importance 
of this point. 

B. Conceptual difference between Wannier and Frenkel excitons 

1. Wannier excitons 

Wannier excitons are constructed on delocalized electrons excited from the valence band to the conduction band. 
These semiconductor bands result from the periodic ionic structure of the semiconductor lattice. A simple way to 
make these bands appearing is to add and substract a one-electron operator 



V ee = J2 7 "ee{r i ), (2.5) 



to the semiconductor hamiltonian H. Yet arbitrary, V ee physically represents a mean electron-electron interaction. 
We will show below the appropriate way to choose it. 

This leads us to rewrite the semiconductor hamiltonian H, given by Eq. (2.1), as 

H = H ( ™ ] +V couU (2.6) 
where V cou i, usually called "semiconductor Coulomb interaction", corresponds to the difference 

Vcoul ^ee ^eei 

The zero order hamiltonian for Wannier excitons Hq W ^ is a sum of one-electron operators. It can thus be written as 

Hq ^ = Hki n ~\~ Ve—ion ~\~ Vee + Vion — ion = ^ ' h\ \ (2-7) 
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where the one-electron operator h\ ' is given by 

AT = ^+E + *«(r,) + i-y io „_ io „ = ^ + „(r,). (2.8) 

This one-electron operator has the lattice periodicity if v ee (r) is chosen with such a periodicity. Besides this require- 
ment, we must also enforce v ee (r) to be such that the resulting interaction v(r) defined in Eq. (2.8) fulfills 



J dr v(r) = 0, (2.9) 

in order for h^ w ^ to be the hamiltonian of a fully neutral system. The simplest choice for v ee (r) is to take it as a 
constant, through the so-called "positive jellium", namely, v ee (r) — N~ Vi on -i on . 

Due to the periodicity of the potential v(r), the eigenstates of hS w ^ are made of delocalized states separated by band 
gaps. The relevant ones for the physics of Wannier excitons belong to the last filled band, called valence band, and 
the first empty band, called conduction band. Close to these band extrema, the eigenstate energies can be written as 
A + fi 2 k 2 /2ra c and h 2 k 2 /2m v , where A is the band gap, while m c and m v are the electron masses for the conductance 
and valence bands dressed by the lattice periodic potential. Note that the effective mass for valence electrons, which 
are close to a maximum is negative; so that resulting valence hole mass, defined as nih = —m v , is positive. 

The eigenstates \vk) of the hamiltonian h^ w ^ with v equal to v or c for valence and conduction states 

{h<W-e vk )\vk)=0 (2.10) 
are then used as a one-electron basis, to rewrite the semiconductor hamiltonian H in the second quantization. 



2. Frenkel excitons 



The situation for Frenkel excitons is totally different: while Wannier excitons are constructed on delocalized valence 
and conduction electron states, the physical picture of the semiconductor excitations giving rise to Frenkel excitons is 
a set of electrons tight to their ions, these electrons switching from the atomic ground state to the atomic first excited 
level. Consequently, the physically relevant one-electron states for Frenkel excitons are not the free (delocalized) 
electrons in a periodic lattice, used in the case of Wannier excitons, but instead the electron localized atomic states 
associated to the various ion sites. 

In order to make these physically relevant atomic states appearing, we are led to perform a grouping of terms in 
the semiconductor hamiltonian H , given in Eq. (2.1), different from the one we have done for Wannier excitons. This 
new grouping of terms is 

H = Hq F ^ + V e - e + Vion-ion. (2-H) 
( F) 

The zero order hamiltonian for Frenkel excitons Hq still is a one-electron operator, but it now contains the electron 



excitons Hq W K It precisely reads 



kinetic contribution plus the electron-ion potential; so that it differs from the zero order hamiltonian for Wannier 

Hff ^ = Hhin + V e -i on = ^2 hi, (2-12) 



where hi is now given by 



„2 



=i 



C. Atomic states 



We can note that, in the one-electron hamiltonian h il enter the interactions of the electron i with all the ions 
n; so that hi differs from a simple atomic hamiltonian. Nevertheless, it is rather clear that the physically relevant 
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one-electron states for Frenkel excitons are going to be these atomic states, i.e., the eigenstates of one electron in 
presence of one ion. Let us introduce them. 

(i) The atomic states \v) for one ion located at R = 0, associated to the Hamiltonian 

2 2 

hatom = £---, (2.14) 

2m r 

are such that 

(h atom -e v )\v)=0, (2.15) 
their wave functions being <p v (r) = (r\u). As the Hamiltonian eigenstates form an orthogonal set, we do have 

(f» = f dript, (r)<p u (r) = <W (2.16) 

(ii) If we still consider one ion but located at R n instead of R = 0, the corresponding atomic Hamiltonian reads 

h(n) = t~- i C \ r (2-17) 
2m |r-R n | v ; 

Due to translational invariance, the eigenstates of read in terms of the atomic hamiltonian states (e v , as 

{hW-e v )\vn) = 0, 

their wave function being such that 

(r\vn) = ipvn{r) = vA v ~ R «)- (2-18) 

For a given ion R n , the bound and extended states of the hamiltonian form a complete basis for one-electron 
states; so that we do have 

(ni/\m) = <W, (2.19a) 

I = ^2\m)(nv\. (2.19b) 

V 

with the sum restricted to the atomic levels v. 

( F) 

(iii) If we now turn to the hamiltonian Hq , given in Eq. (2.13), we see that it differs from a bare sum of atomic 
hamiltonians, since each electron feels the interaction of all the other ions. This is a real difficulty: unlike Wannier 
excitons, in which the eigenstates of Hq W ^ can be obtained exactly, the diagonalization of the one-body part of Hq F ^ 
for Frenkel excitons can only be approximated due to this multiple ion interaction. As shown below, this will force 
us to make assumption on the atomic wave functions extensions, the Frenkel exciton picture being appropriate when 
the tight-binding approximation is valid. 

(F) 

Due to the sum over all ions contained in the one-body hamiltonian Hq given in Eq. (2.13), it is, on the one 
hand, clear that the states \vn), for a fixed n (which form a complete set for one-electron states due to Eq. (2.19)) 
cannot be a physically relevant basis to describe Frenkel excitons made of excitations on all possible ion sites. On the 
other hand, it is also clear that if we now leave n running over all the ion positions, the states \vn) for all v and all n 
form an overcomplete set - the states \vn) for one particular n forming a complete set already. As a bare consequence 
of this overcompleteness, the states \un) with different n are not exactly orthogonal. 

However, when the atomic states v and v' are highly localized compared to the ion-ion distance, the wave function 
overlap between different atoms is essentially zero 
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<pt,(T - R n >K(r - R n ) ~ for n' £ n, (2.20) 
whatever r is. Consequently, the scalar product of such atomic states reduces to 

(n'u'\uri) = J drtp*,(r - R„/)^(r - R„) ~ <W<W (2.21) 

due to Eq. (2.16): These states \vn) are thus quasi-orthogonal. 

(iv) In spite of these difficulties, linked to the overcompleteness of the states \vn) with n running over all ion sites, 
the \vn) states have to play a role in the Frenkel exciton physics. This is why we are going to use them in the second 
quantization description of Frenkel excitons. 

Before going further, let us add some comments on using these states \un) as a basis for second quantization. It 
is clear that Eqs. (2.20) is not valid for atomic extended states, nor even for the highest bound levels. The states 
\vn) can however be seen as a nice basis for one-electron states if the problem at hand relies on the highly localized 
(lowest) atomic states, for which the overlaps between atomic wave functions for different ions is negligible. This is 
actually the case for conventional Frenkel excitons, in which the v's of interest reduce to the ground state [y = 0) and 
the first excited state {y = 1). 

A somewhat cleaner way to present the use of these overcomplete states \vn) is to say that we can always add 
to the two sets of atomic states of physical interest, namely, \v = 0, n) and \v = l,n), other states constructed "in 
an appropriate way" , in order to form a complete orthogonal basis when added to the two states \v = (0, 1), n). In 
problems physically controlled by the two lowest atomic states, these additional "appropriate states" are not going 
to play a role in the final results. This is why it is far simpler not to consider them at all and to stay with the full 
overcomplete set of states \un) for all n and all v, the i/s different from (0, 1) playing no role in the end. As a direct 
consequence, in the following, the sums over v will have to be considered as sums over v = (0, 1). 



III. THE FRENKEL EXCITON HAMILTONIAN IN THE SECOND QUANTIZATION 



A. One-electron creation operators relevant to Frenkel excitons 



Let a\, n be the creation operator for the atomic state \vn), 

\vn) = at n \v) . 



(3.1) 



Since the state \vn) expands on the plane wave basis |k) as \vn) = ^ k |k) (k\vn) , the atomic state creation operator 
a\, n reads in terms of the free electron creation operator at as 



(3.2) 



By using the fact that the electron operators anticommute, 



— ^kk'j h is easy to show, using the above 



relation, that the anticommutator for atomic state destruction operators is exactly zero 
while for the lowest highly localized states v = (0, 1), we do have, due to Eq. (2.21), 



(3.3) 




B. One-body part of the Frenkel exciton hamiltonian 

Let us now use these operators to rewrite the part Hq F ^ of the semiconductor hamiltonian given by Eq. (2.12). 
Since is a sum of one-body operators, it can be written in terms of the creation operators for the one-electron 
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states \vn) relevant to Frenkel excitons, as 



H o F) = E £ <V n '' ^)4'n'<W 



(3.4) 



v n 



According to the second quantization procedure, the prefactor s(v'n' vn) is given by 



e{v'n',vn) = J dr ^,„,(r) ^ - £ 
Due to Eq. (2.18), this prefactor also reads 

/2 
dr <p*, n , (r) e„ - \ r 



<A,„(r). 



¥>™(r). 



(3.5) 



(3.6) 



So that, for highly localized states, as the ones of physical interest, it reduces, due to Eq. (2.20), to 

e{v'n' , vn) ~ 5 nn < [e y 5 vv > + v{v' , v)] , (3.7) 

where v(y\ v) comes from the interactions with all the other ions. Due to the translational invariancc of the hamil- 
tonian leading to Eq. (2.18), v{v' ', v) can be rewritten as 



(,» = £ [dr rt,(r) £ ■ rR e2 R ^ W = (-'I E i^R> 



(3.8 



where the vectors R in the sum correspond to all possible distances between ions. The one-body part of the semicon- 
ductor hamiltonian appropriate to Frenkel excitons ends by reading as 



H F) = E £ « a ln a ™ + E V ( V ' ' V ) a l'n a » 
v,n n,v'^v 



(3.9) 



where e v = e v + v(v, v). The second term in describes the fact that the electron of a given site can change its 

atomic level from v to v', while staying on the same site, due to its interaction with the ions of the other sites, as 
seen from the definition of v (i/, v), given by Eqs. (3.8). However, for the states v — (0, 1) highly localized compared 
to the ion-ion distance, these v{y' ', v) scatterings are extremely small. 



C. Electron-electron interaction 

Let us now turn to the two-body operator V eei defined in Eq. (2.4). The standard second quantization procedure 
leads to write it on the atomic basis \vn) as 

where the prefactor is given by 

For v and v' equal to (0, 1), this prefactor is nonzero for n[ — ri\ and n' 2 = n 2 only, due to Eq. (2.20). If we then 
use the translational invariance of atomic wave functions, namely, Eq. (2.18), it is easy to see that V ee can be written 
as 

Vee = \ V R ni -R„ 2 [y x H ) a l' in A' 2 n 2 a ^ 2 a Vini (3.12) 

niv[vi 
ri2 v' 2 vi 

where the electron-electron scattering depends on the distance R between ions through 

Vn ( 4 H ) = V- R ( U \l\)=j Mr^t, (r lK , (r 2 ) ^ _ ^ + R| ^ (r^ ( ri ). (3.13) 
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D. Restricted hamiltonian for Frenkel excitons 

The atomic states [y = 0, n) and [y = 1, n) are the equivalents of the valence and the conduction band states for 
Wannier exciton. By noting that the energy of the atomic ground state e„ = o is very different from the one of the first 
excited level £„=i, we are led to think that the physically relevant part of the hamiltonian corresponds to processes 
in which the number of electrons in the v — level and the number of electrons in the v = 1 level are separately 
conserved. This is equivalent to processes which keep the number of conduction electrons and the number of valence 
electrons fixed, for the case of Wannier excitons. This leads us to drop terms with v ^ v' in Hq F \ The latter then 
reduces to the first sum of Eq. (3.9), namely 

H f ] - H o =£o^2 a o™ a o« + £i ^2 a L°in- ( 3 - 14 ) 

n n 

If we now turn to the Coulomb interaction V ee , given in Eq. (3.12), and we also keep terms which conserve the 
number of v = and v = 1 electrons separately, we see that they are of two kinds, (i) V ee contains intraatomic 
processes, in which the electron of a given site stays on the same level. These processes correspond to terms like 
al ni al n2 a vn2 a vni with v equal or 1. V ee also contains terms like ao ni a\ n2 ain 2 a On 1 with and 1 possibly exchanged 
- which makes such a term appearing with a factor of 2. (ii) In addition, V ee contains interatomic processes, in which 
the electron of one site jumps from v — to v — 1, while the electron of another site goes from v = 1 to v = 0, 
namely, terms like al ni a\ n2 a 0ri2 ai ni , with and 1 possibly exchanged - which makes this term also appearing with a 
factor of 2. 

So that V ee for Frenkel exciton ultimately reduces to V ee — vjP + Ve? + + V^e with 



v ee 



y( 2 ) = 

v ee 



v ee 



v ee 

Note that the last two terms differ from zero for n\ — n 2 , while the first two terms are equal to zero for n\ = n 2 ; 
this is why we have excluded n\ — n<i from the first two sums. 
These four terms are shown in Fig. 1. 

IV. ELECTRON-HOLE HAMILTONIAN FOR FRENKEL EXCITONS 
A. Electron and hole creation operators 

As for Wannier excitons, it is appropriate to introduce the concept of hole. This will allow us to start with a |0) 
state, in which the electrons of all sites are in the atomic ground state, v = 0, and to speak in terms of excitations with 
respect to this ground state, i.e., in terms of the small number of sites in which the electron is no more in the ground 
state, this number being one for one exciton, two for two excitons and so on... Such an elementary electron-hole 
excitation is shown in Fig. 2. Due to the electron-hole attraction and the cost in electrostatic energy induced by 
the electron and hole separation, we expect the lowest energy excited states to correspond to n\ = n 2 . As a direct 
consequence, the Frenkel excitons are going to be made from electron-hole pairs on the same site. Let us now recover 
this obvious result. 



1 t/ f° °\ t t 

9 2^ K ^i-«»2 1 ) a 0n 1 a 0n 2 a 0n2 a 0n 1 , 



1 



1 1 



^2 V Rm-Rn 2 ( ) a 0n 1 a ln 2 a ln 2 a n 1 



1/ f 1 \ t t 

VR ni -Rn 2 I 1 q I a l„ 1 a 0«2 a l«2 a 0n 1 - 



(3.15) 



9 



In the absence of spin degrees of freedom, the electron and hole creation operators are simply linked to the v = 
and v = 1 atomic state level creation operators through 

An = 4 

a 0n = b\. (4.1) 

By using the anticommutation relations for electrons in atomic states, given by Eq. (3.3), it is straightforward to 
show that [a n ',a n } + — 0, while 

[a n >,a n ] + = a ln ',a{ n = (ln'\ln) ~ S nn > (4.2) 

for highly localized atomic states compared to the interatomic distance, which makes Eq. (2.21) valid. In the same 
way, [b n ,,b n ] + = 0, while 

[bn',bl] + = [aj n ,,a „] = (0n'|0n) ~ <W (4.3) 

If we now turn to the anticommutator between electron and hole operators, we find that 

[ a n',bl] + = [a\ n >,a 0n ] + = 0, 

[a n ',b n ] + = [ai n /,aj n + = (ln'\0n) ~ 0, (4.4) 
for highly localized atomic states, this last anticommutator being exactly equal to zero for n — n', due to Eq. (2.19a). 

B. One-body operator Ho in terms of electrons and holes 

In order to rewrite the part Hq of the Frenkel exciton hamiltonian given in Eq. (3.14) in terms of electron and hole 
operators, we first note that aJ n aon = 1 — a 0n% n = 1 — b^bn, so that 

H = N s e + ^(eiolan - e ^&n)- (4.5) 

n 

Due to other contributions in ajja„ and 6jj6„ coming from the electron-electron interaction, the electron energy and 
the hole energy are going to differ from E\ and (— e ), as now shown. 

C. Electron-electron interaction in terms of electrons and holes 

We now turn to the electron-electron interactions given in Eq. (3.15). Since a\ n = a n , the second term of this 
equation, shown in Fig. 3a, readily gives the Coulomb repulsion between two electrons as 

K ( e 2) = V ee = \ ]T V R ^-R n2 (\\\ 4 2 fl„ 2 O ni . (4.6) 

In order to rewrite the third term of Eq. (3.15), we first note that aj ni a| n2 a 1 „ 2 a „i = (— o\n 2 a 0ni)(~ a a ni a in 2 )7 
while aJ ni ao, u = 1 — ao ni a\ n ^\ so that this third term gives two contributions, 

k^k^+E^E^-^ (J I)- ( 4 - 7 ) 

«2 "1 

(3) 

The second term of Vee , which comes from the Coulomb interaction between one electron in an atomic excited level 



of Vee , given by 



and all the atomic ground states, is going to dress the bare electron energy e\, appearing in Eq. (4.5). The first term 

Vef r) = E *V-«« 2 (II) biA 2 a n2 b ni , (4.8) 
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describes a direct electron- hole attraction, the electron and the hole staying in their sites (see Fig (3c)). 

In the same way, the fourth term of Eq. (3.15) gives two contributions since a \n 1 a on 2 ain 2 a On 1 = 
a L 1 a On 2 (~ a o«i a in2) wnile - a on 2 a Om = -^112 + a o ni al n2 ; so that it reads 

= V (e xch) _ Vr= J 1 \j2ala n . (4.9) 

The second term of vj^ is also going to dress the electron energy e\, while the first term, given by 

Vt ch) = E V ^-Rn 2 ( J J ) <blb n2 a n2 (4.10) 



and shown in Fig. (3d), corresponds to the destruction of one electron- hole pair on the site n<i and to its recreation 



on the site n\. Let us stress that, while the direct electron-hole Coulomb interaction V^ lr ^ in Eq. (4.8) corresponds 



to an attraction, the exchange electron-hole Coulomb interaction V^ xch ^ is repulsive. 

We now turn to the first term of Eq. (3.15) between two ground state electrons v = 0. We first note that, due to 
Eqs. (3.3-4), 

a 0n 1 a ln 2 a 0n2 a 0n 1 = 1 - a0ru a ni - a 0n 2 a 0„ 2 + a 0n 2 «0ru aj m a 0n 2 ■ (4-11) 

So that, by grouping the two terms with a minus sign, VeP generates three contributions 



The last term of vjP, shown in Fig. (3b), corresponds to an hole- hole repulsion. It precisely reads 



V hh = \ E V^-Rn, J! blbt 2 b n2 b ni . (4.13) 

rii#n 2 



The first term of vjP is a bare constant which describes all Coulomb interactions between ground state atomic levels. 
It produces a band gap renormalization. The second term of VeP comes from the interaction between one particular 
ground state electron, in the site n, and the other ground state electrons. This term has to appear when the site n 
is empty, i.e., occupied by a hole, in order to compensate for the electron-electron interaction already included in the 
constant term of VeP- This second term is going to dress the atomic ground state energy Sq, when speaking in terms 
of holes, as fully reasonable, since all interactions between atomic ground state electrons are by construction forgotten 
when we turn to electrons and holes. These interactions actually appear through the renormalization of the atomic 
ground state and excited state energies, the electron and the hole being more subtle objects that just one electron in 
the atomic excited state and one electron absence in the atomic ground state. 

D. Electron- hole hamiltonian 

If we now collect all these terms, we end by writing the part of the semiconductor hamiltonian, appropriate to 
Frcnkcl excitons H + vjP + v£ ] + v£ ] + V e ( e 4) + V ion - lon as 

=A + H eh + V mtra + V mter . (4.14) 
(i) A is a constant which contains contributions from the atomic level ground states only. It precisely reads 

A = N s s + N s v(0,0) + l E *V-*» a (S o)+^ E | Rw | -N s (e + s^ l) ) (4.15) 

ni#n 2 ^ / ni^n 2 ni 
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where N s is the number of ion sites. By using Eqs. (3.8) and (3.13), the Coulomb contribution to this band gap 
renormalization, given by the bracket of the above expression, can be rewritten as 



4 COU ° = E [drdr'\Mr)\ 2 \Mr') 



1 



|r-R| 2|r-r'-R| 



If! 

2 R 



(4.16) 



where the R's are the possible distances between two ions. Note that the last term, e 2 /R, which comes from the 
ion-ion interaction, and which makes the system at hand neutral, allows for the convergence of in the large 

sample limit. 

(ii) The second term of Eq. (4.14) is a one-body operator which can be written as 



H e h = e e £ a n a n + e h £ b\b n . 



(4.17) 



It describes the electron and hole kinetic energies. These energies, given by 

R V 







V R=0 



1 

1 



(4.18) 



e h = £ + v{0, 0) + V V R I " " ) (4.19) 







^ K \ 

differ from the atomic bare ground and excited state energies — £o and e\ due to Coulomb interactions with all 
the atomic ground states. These contributions have to appear, when we speak in terms of holes since all Coulomb 
interactions among these atomic ground state levels are then forgotten, by construction. 

(iii) The third term Vintra of Eq. (4.14) corresponds to V^ zr ^ + V^ xch ^ taken for n\ = n 2 - It precisely reads 

V intr a = S'y] a* 6 n a n , (4.20) 



where (—5) is given by 



= » ( I I )+Vr.=o( J I ) (4.21, 



By using Coulomb couplings, given in Eq. (3.13), we see that this quantity also reads 

S = J dndra [^(r 2 )^(n) - ^(r 2 )^(n)] y*(r 2 )yS(n). (4.22) 

This shows that 5 is a positive constant, since {v\v) — 1, while (0|1) = 0. This energy 5 physically corresponds to the 
energy decrease when the site n is occupied by an electron- hole pair, i.e., when the site n is neutral. This is going to 
make the potential Vintra responsible for the fact that excitons are made from linear combinations of electrons and 
holes located on the same site. 

(iv) The last term of Eq. (4.14), Vinter, is made of all Coulomb interactions between sites. It contains the electron- 
electron and hole-hole contributions V ee and Vhh, given by Eqs. (4.6) and (4.13), which are interactions between 
sites by construction, since a given site cannot accommodate two electrons or two holes due to the Pauli exclusion 
principle. It also contains the part of the direct electron-hole potential V^ lr \ taken for n\ ^ n 2 . Using Eq. (4.8), 
this direct electron-hole exchange interaction between sites precisely reads 

V e f r) = ~ £ V Rni . Rn2 (J J) blala n2 b ni . (4.23) 

It finally contains the part of the electron-hole exchange potential V^ xch ^ given in Eq. (4.10), taken for m ^ n 2 - 
This part has a very special role since it allows the excitation transfer from one site to the other. Let us isolate this 
transfer term from the other Coulomb terms and call it as Vt ra ns 

V trans = £ V R ni ~Rn 2 f J J J «UUni«nr ( 4 - 24 ) 
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All this leads us to write Vinter in Eq. (4.14) as 

Vinter — Vtrans ~\~ Vcoul j 

Vcoui = V ee + V hh + V e f r \ (4.25) 
The four contributions of this Vi nte r operator are shown in Fig. (4). 

E. Discussion 

The expression of the semiconductor hamiltonian appropriate to Frenkel exciton in terms of electrons and 

holes, given in Eq. (4.14), allows an easy comparison between highly localized states leading to Frenkel excitons, and 
extended states, leading to Wannier excitons. We first see that the electron and the hole energies for Frenkel excitons 
appearing in H e h are constant, while they depend on k for Wannier exciton: the electrons and holes for Wannier 
excitons, which belong to the conduction and valence bands, are delocalized over the whole sample, so that their 
energies must depend on momentum. In addition, the " electon-hole exchange", i.e. the possibility for one electron- 
hole pair to recombine while another pair is created, plays essentially no role for Wannier excitons: it is just responsible 
for a small splitting between Wannier excitons, when the spin degrees of freedom are included. On the opposite, this 
"electron- hole exchange" is crucial in the case of highly localized states, as it is the only process allowing an excitation 
transfer between sites: this makes the operator Vtrans entirely responsible for the Frenkel exciton formation, as we 
now show. 



V. LOWEST EXCITED STATES IN THE ABSENCE OF INTERACTION BETWEEN SITES 

Let us first forget the interactions between the sites 

H-pair - 

if we drop the irrelevant band gap renormalization A. 



. The hamiltonian H^ F > , given in Eq. (4.14), then reduces to 

= H e h + Vintra , (5-1) 



A. Ground state and lowest excited states 



The ground state of H pa i r has zero electron-hole pair. Let us call it |0) and take its energy as 0. 

If we now consider the one electron-hole pair state a^b^, |0) with an electron on site n and a hole on site n', we see 
that its energy is e e + Eh for n ^ n' while it is e e + Sh — 5 for n = n' . Since <5 is positive, the lowest excited states of 
H P air thus have one electron-hole pair on the same site. They reads 

(Hpair -^pair) \Rn) 0? 

|i?„) = 4^|0)= J Bt| ) (5.2) 
with Ep a i r = e e + Eh — S. These states form a ^-degenerate subspace, since n can run from 1 to N s . 



B. Commutation rules 



Using the anticommutation rules for electrons and holes given in Eqs. (4.2-4), it is easy to show that the electron-hole 
pair operators B^ behave as bosons with respect to the destruction operators since their commutator reads 

[B n ,B n ,]_=0, (5.3) 

while they are composite bosons only since the other commutator is such that 

[B n ,,Bl]_=5 nn ,-D nn ,, (5.4) 

the deviation-from-boson operator for electron-hole pairs being equal to 

Dnn> = <W {ala n + b^bn). (5.5) 

As standard for deviation-from-boson operator, D nn i gives when acting on the electron-hole pair vacuum. 
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VI. FRENKEL EXCITONS 

If we now keep the coupling between sites Vi nte r — V tra ns + V cou i, denned in Eq. (4.25), in the Frenkel exciton 
hamiltonian H^ F \ we induce non-diagonal contributions between different sites. They are going to split the degenerate 
subspace \R n ). The Frenkel excitons result from the diagonalization of the hamiltonian in this \R n ) degenerate 

subspace. 

A. Derivation of the Frenkel excitons 

Since the states \R n ) have one electron-hole pair only, the electron-electron and hole-hole parts V ee and Vhh of 
Vcoui, defined in Eq. (4.25), give zero when acting on \R n )- The direct electron-hole interaction V^ tr \ given in Eq. 
(4.23), also gives zero since b n2 a ni \R n ) = for ri\ ^ n 2 . Consequently V cou i \R n ) = 0; so that the only part of Vi nter , 
which plays a role in the diagonalization of H ( F > in the \R n ) subspace is the electron- hole exchange term Vtrans- Since 
b„ 2 a n2 al bl |0) = S ri2n |0), we readily find 

Vtrans \R n ) = E V Rn ^ Rn ( J J ] \R ni ) . (6.1) 

This shows that, if we drop the irrelevant constant A, the Frenkel exciton hamiltonian H ( F ) acting on \R n ) reduces 



to 



tr(0) tt i T/ 

11 X — 1A pair t vtrans? 



H^x \Rn) — E palr \R n ) + E VR ni -R n | j j| j \Rni) ■ 



(6.2) 



Let us now show that the following linear combinations of \R n ) 

\X Q ) = -!=Y,e^\R n ), (6.3) 



v s n=l 

known as Frenkel excitons, are the exact eigenstatcs of the hamiltonian To prove it, we first consider V tr ans 

acting on \Xq) 



Vtrans \X Q ) = -j= £ ^ ^ £ (\ \\ 



(6.4) 



which also reads 



Vtrans , v • .. , . 

n' n^n' 

In order to calculate the last sum, we rewrite it as 



\X Q ) = 1 £ \Rn>) E ^ Q R "^„^„, ( J I ) ■ (6.5) 



(6.6) 



Due to the invariance of the system, the sum in RHS cannot depend on n'\ so that the above equation leads to 



E ^V Rn ,. Rn (; J)= e**"* E e " iQR ^ ( ? 



1 




(6.7) 



where the sum is taken over all possible distances R between ions. When inserted into Eq. (6.5), this readily leads to 

Vrans \X Q ) = \X Q ) E e^^R (\\\ (6-8) 
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So that we end with 

Hf \X Q ) = E Q \X Q ) , (6.9) 

where the eigenenergy is given by 

E Q = E pair + ]T e~ iQIl V R (H). (6.10) 

The above equation shows that the splitting of the N s degenerate states \R n ) into N s states \Xq) is only due to 
the electron-hole exchange V tr ans between different sites: V tr ans is the only part of the hamiltonian allowing for an 
excitation transfer from site to site, as necessary for the derealization of the Frcnkcl excitons. 



B. Frenkel exciton creation operator 



Eq. (6.3) leads to write the Frenkel exciton creation operator Bq defined as \Xq) — Bq |0) in terms of the creation 
operators B^ — a^b^ for one electron-hole pair on site n as 



N 

Q -fN s 

n=l 



= Ee iQR "5t. (6 . n) 



In the same way as free electron-hole pairs can be written in terms of Wannier excitons, it is possible to write 
electron-hole pair localized on site n in terms of Frenkel excitons. This is barely done by noting that 

=E^ qr " b q = w E it E e * Q(R "'- Rn) = B i (6-12) 



V s Q 6 n'=l Q 

since the sum over Q is equal to for n' ^ n and to N s for n' = n. 

These Frenkel excitons are expected to be composite bosons. This is easily seen from their commutation rules. By 
using Eq. (5.3) for electron-hole pairs on site n, we readily find 

[B Q ,B QI ]_=0. (6.13) 

If we now turn to the other commutator, it reads 

- = w E E e " iQ ' Rn ' eiQR " [^iBi]-- (6-i4) 

n' — l n—1 

We then use Eq. (5.4) for the commutator \B n i,B^\_. The term in 5 nn i gives N~ 1 ^ n e l ( Ci ~ < ~ 1 ) Rn , which is equal 
to for Q ^ Q' and 1 for Q = Q'; so that we end with 



B Q' ) B Q 



Sq>q - Dq'q (6.15) 



where the deviation- from-boson operator -Dq'q for Frenkel excitons, which comes from the deviation- from-boson 
operator for electron-hole pairs D n i n appearing in \B n i,B^\_ reads as 

Dq'q = w E eJ(Q_Q ) R "(4«« + bib n ). (6.16) 

s n—l 

We can note that, as for Wannier excitons, this deviation-from-boson operator gives when acting on the pair 
vacuum state 

£>q' Q |0)=0. (6.17) 

This leads us to conclude that, in order to describe the interactions between Frenkel excitons properly, it is necessary 
to follow a path similar to the one we have used for Wannier excitons, namely, to define the Pauli scatterings of two 
Frenkel excitons for carrier exchanges without Coulomb interaction and the Coulomb scatterings of two Frenkel 
excitons for carrier interactions without carrier exchange. As for Wannier exciton, the composite nature of the 
Frenkel excitons makes impossible a clean description of the interactions between two excitons as a potential, the 
only well-defined quantity again being the "creation potential" of the Q exciton. The calculation of these Pauli and 
Coulomb scatterings for Frenkel excitons, necessary to handle their many-body physics properly, will be done in the 
forthcoming publication. 
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C. Interacting Frenkel exciton hamiltonian 



If we have more than one electron-hole pair, the Coulomb part V cou i of the Frenkel exciton hamiltonian, given in 
Eq. (4.23), is going to play a role. This leads us to rewrite as = H x with 



where H x , as given in Eq. (6.2), corresponds to the pair hamiltonian H pair plus the part of the Coulomb interaction 
Vtrans allowing for the excitation transfer. 

The remaining part V cou i, given in Eq. (4.25), corresponds to all direct Coulomb interactions between two electrons, 
two holes, and one electron and one hole in different sites. This operator is going to generate all many-body effects 
between excitons induced by Coulomb interactions. In addition to them, as for Wannier excitons, also exist many-body 
effects induced by Pauli exclusion through the fact that Frenkel excitons are not elementary bosons. The scatterings 
associated to V cou i and to the deviation-from-boson operators will be calculated in a forthcoming publication. 



In this paper, we have derived the creation operator for Frenkel exciton starting from the microscopic hamiltonian 
for free electrons in a periodic lattice. Let us summarize the main steps of the derivation. 

1) We first isolate one ion located on site n and we introduce the atomic states \vn), eigenstates for this particular 
ion. They form a complete basis for one-electron states. 

2) If we let n running over all sites, the states \vn) form an overcomplete set. However, if the states of physical 
interest are the two lowest atomic levels, v — and v = 1, the states \vn) with v = (0, 1) are essentially orthogonal 
in the tight binding limit, i.e., when the overlap of the v = (0, 1) wave functions on different sites is negligible. This 
allows us to use the \v = (0, 1) n) as a one-electron basis to describe Frenkel exciton in second quantization. 

3) The large energy difference between the atomic states v = (0, 1) leads us, in the electron-electron interaction 
written in terms of the creation operators a) vn for these \vri) states, to only keep the terms which conserve the number 
of electrons in state v = and in state v = 1 separately. 

4) When written in terms of electron-hole pairs with a* = a[ n and b\ L — a 0n , the electron-electron interaction 
generates a constant term, which renormalizes the band gap. It also generates one-body contributions in 6^6 n and 
ajjdn which dress the (y = 0) and (y = 1) atomic levels, the electron and hole energies differing from the atomic 
energies due to the Coulomb interaction with a kind of jellium having one electron in the ground state of each ion 
site. 

5) Finally, the electron-electron interaction also generates an intersite contribution V%nter given in Eq. (4.22a) and 
an intrasite contribution Vintra which insures local neutrality; so that the lowest excited states have one electron and 
one hole on the same site. 

6) The corresponding degenerate subspace a^b^ |0) = |0) with n running over all sites, is split by the part Vtrans 
given by Eq. (4.21) of the intersite contribution Vmter which allows to transfer the excitation from one site to the 
other. 

7) The resulting eigenstates correspond to a set of dclocalized excitations, known as Frenkel excitons. Their creation 
operators read 



where N s is the number of ion sites, these ions being located at R„ on a periodic lattice. 

8) These Frenkel excitons are composite bosons. Their many-body effects thus have to be handled along a procedure 
similar to the one we have used for Wannier excitons. We are going to show that they predominantly interact through 
the Pauli exclusion principle between their electron-hole components which make them differing from elementary 
bosons and which produces "Pauli scatterings" to describe carrier exchanges without carrier interaction. Frenkel 
excitons also interact through the V cou i part of the intersite interaction given in Eq. (4.25), which contains electron- 
electron, hole-hole and electron-hole direct processes, i.e., processes in which the carriers stay on their site. 

The composite-boson many-body theory appropriate to Frenkel excitons will be presented in a forthcoming publi- 
cation, the present work providing the necessary tools to build this theory on solid grounds. 



H x = H ( p + V c 



coul 1 



(6.18) 



VII. CONCLUSION 
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VIII. FIGURE CAPTIONS 



Fig 1. Electron-electron potentials, given in Eq. (3.15). 

(a) Terms in aQ„ i aQ„ 2 aon 2 a On 1 , in which the two electrons stay in their atomic ground state v = 0. 

(b) Terms in a\ ni a\ n2 ai n2 a\ ni , in which the two electrons stay in their atomic excited state v = 1. 

(c) Terms in al ni a\ n2 ai n2 a 0ni , in which one electron stays in the ground state, while the another one stays in the 
excited state. 

(d) Terms in a\ ni al n2 ai n2 av ni , in which the electron on site n\ is excited from the ground state v = to the excited 
state v = 1, while the electron on site n-i returns to its ground state. 

In these four processes, the number of electrons in the ground state and in the excited state are separately conserved. 
Fig 2. Excitation of an electron from the ground state on site n\ to the excited state on site n-i'- this corresponds 
to the creation of an electron-hole pair on sites (ri2, n\). 

Fig. 3. (a) Electron-electron interaction V ee , given in Eq. (4.6). 

(b) Hole-hole interaction Vhh, given in Eq. (4.13). 

(c) Direct electron- hole interaction V^ lr \ given in Eq. (4.8). 

(d) Exchange electron- hole interaction V^ xch \ given in Eq. (4.10). 

Fig. 4. Interaction between sites, described by Vi nter , given in Eq. (4.25). 

(a) The part V tra ns, given in Eq. (4.24), describes the destruction of an electron- hole pair on site n\ and its 
recreation on site n 2 . 

(b) Direct interactions between two electrons, corresponding to V ee , given in Eq. (4.6), (c) between two holes, 
corresponding to Vhh, given in Eq. (4.13), and (d) between one electron and one hole, corresponding to , given 
inEq. (4.23). 



